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We derive WKB expressions for glueball masses of various finite temperature super- 
gravity models. The results are very close to recent numerical computations. We argue 
that the spectra has some universality that depends only on the dimension of the AdS 
space and the singularity structure of the horizon. This explains the stability of the ++ 
glueball mass ratios between various models. We also consider the recently proposed non- 
supersymmetric model arising from the type string. In the supergravity limit of this 
model, the heavy quark potential has an effective coupling with 1 / (log u) behavior in the 
UV. Unfortunately, the supergravity solution implies that the heavy quark potential is still 
coulombic in the infrared, with an effective coupling of order 1. We also argue that the 
type supergravity background solution does not have normalizable glueball solutions. 
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1. Introduction 



One of the many interesting developments to arise out of Maldacena's conjecture 
[|2]||] is the ability to study nonsupersymmetric large iV gauge theories at strong coupling 
[§-|J. One studies a d dimensional euclidean gauge theory at finite temperature, which is 
equivalent to a theory with d — 1 noncompact directions and a Euclidean time compactified 
on a circle of circumference (3. As was pointed out by Witten f|] the Maldacena conjecture 
relates wave equations in an AdS blackhole background to two point correlation functions 
of a finite temperature Yang-Mills gauge theory. 

Using this conjecture Witten argued that the dilaton wave equation in this background 
implies a discrete glueball spectrum with a finite gap. This spectrum was studied more 
closely in and also in fl(],|ll[ where comparisons were made between the supergravity 
results and lattice gauge theory results. However, the strong coupling behavior of QCD 
is highly nonuniversal, so there really is no reason to expect much similarity between 
the lattice results and the supergravity results, beyond the fact that the spectra for both 
theories is discrete with a finite gap. 

However, one might hope to find some universality within different supergravity mod- 
els. In particular, other supergravity models were recently studied that correspond to 
finite temperature QCD with its R symmetry group broken [p~2] , p~3f1 . This has the nice fea- 
ture of getting rid of some of the unwanted Kaluza-Klein states. It was noted in [13[ that 
there seemed to be some universality in the mass ratios of J PC = ++ glueball states for 
the different supergravity models. One of the purposes of this paper is to explain this 
universality by finding WKB approximations for the glueball spectra. As it turns out, 
the leading order term depends on the particular supergravity theory being considered. 
However, the subleading term has universal behavior, depending only on the dimension of 
the AdS space and the horizon singularity. One could then speculate that if "real" QCD 
has a supergravity dual that is asymptotically AdS n , then it too will have some universal 
behavior, that depends only on n and the singularity structure at a horizon. 

In all of the above models, supersymmetry was broken by turning on a temperature. 
Therefore, supersymmetry is restored in the ultraviolet and hence the supergravity solu- 
tions do not exhibit asymptotic freedom. Recently, a nonsupersymmetric gauge model 
was proposed that arises from D3 branes in the nonsupersymmetric type theory |T^ , |T5| . 
Since supersymmetry is never restored, one should expect to see running of the coupling 
in the UV. In principle, one could also apply the WKB analyis to this model and derive 
its glueball spectrum, up to an overall scale. 
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The type theory has a tachyon in the bulk that presumably gets an expectation 
value. The authors in [|15| derived a supergravity action for this model. We will find 



asymptotic solutions to the equations of motion coming from this supergravity action. 
The results found here should be taken with a grain of salt since the background metric 
has curvature that is either greater than or roughly equal to the string scale in both the UV 
and IR. Nonetheless, one still hopes that the supergravity results are qualitatively correct, 
as in the case for the entropy of M = 4 super Yang-Mills at finite temperature, where 
the supergravity result differed from the perturbative Yang-Mills result by a factor of 3/4 
||16|| . In fact, in the type case we do find a running of the effective coupling in the UV. 
However, in the infrared we do not find confinement nor do we find normalizable glueball 
solutions. This suggests that one must consider the full a-model to see such behavior. 

In section 2 we derive WKB expressions for the masses of ++ in the finite temper- 
ature models described in [Q. We compare these results to the recent numerical results 
and we find excellent agreement. We also describe the six dimensional model with five 
uncompactified directions and argue that the glueball spectrum has a finite gap with a 
continuous spectrum above the gap. In section 3 we derive WKB expressions for more 
general supergravity models, and we show that these expressions have a general form that 
depends on the dimension of the asymptotic AdS space and the singularity at the horizon. 
This explains the recently noted stability of the spectrum for the class of models discussed 
in [[131 . Using this analysis we also find WKB expressions for the 0~~ and ^ glueballs 
in 3 and 4 dimensions respectively. We again find good agreement with the numerical re- 
sults. In section 4 we describe our findings for the type model. We derive expressions for 
the metric and the coupling in the ultraviolet and infrared and use this to find the heavy 
quark potential in these two limits. We also argue that the IR behavior of the metric and 
coupling does not allow for normalizable glueball solutions. In section 5 we present our 
conclusions. 



2. WKB Masses for 0++ Glueballs 

In this section we compute masses using a WKB approximation for the ++ glueballs 
in the supergravity models in ||. In the next section we will consider more general cases 
as well as the WKB solutions for the other glueballs recently discussed in the literature. 
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The ++ glueball spectrum is governed by the dilaton wave equation in the appropriate 
background. Consider first the background arising fron N stationary Dp branes at finite 
temperature T. The metric in the near horizon limit is 



ds* 



a 



VgN 



1 1 dt 2 + dx 2 



+ 



U 

W 



u(7-p)/2 y 1 _ (ih^y-p 



+ u z dn 



2 jn2 



and with a dilaton background 



g 



jj7-p \ (P-3)/4 



gN 



The temperature T is related to Ut and the coupling by 



T = 



7-p Ut p)/2 



47T 



The dilaton equation of motion is 

d li e- 2 + y /gg' u 'd v <f> = 0. 
Assuming that (f) is of the form cf> = e lk ' x p{u) with k 2 = — M 2 , ( p.4f) reduces to 

du (u 7 - p - U 7 T ~ p ^j Ud u( f> + M 2 gNU(p = 0. 
Defining a new variable x = U 2 and rescaling, ( |2.5| ) reduces further to 



where n = 2/(5 — p) and 



\4TrnT J 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2-5) 



(2.6) 



(2.7) 



The differential equation in (|2.6|) has singularities at x = 0, oo and at all 1/(1 +n) roots 
of unity, so solutions to this equation are unknown for finite n. In order to do the WKB 
approximation, we define a new function ip = * / 2 +i7»_ 4>i an d we change variables to 



x = 1 + e y . The equation in (|2.6|) now takes the form 

( f 



+ \ -e y - -— + - — 



(2- 
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where the primes denote derivatives with respect to y and / = (1 — e y )((e y + — 1). 
For large negative y, the term in front of if) in (|2.8|) is approximately 

For large positive y the asymptotic behavior for this term is 

M^Kf)'-^-*^ ' >>a (2 - 10) 

Thus, for A sufficiently large there will be two turning points at y = — oo and y = yo, where 

y « nlog(4nA/(n + l) 2 ). (2.11) 
Hence, the WKB approximation for this curve gives 



/yo \ If" 1 / f \ 

dydjev - -ij + - \L-J , m > 0. (2.12) 

To leading order in M we may approximate the WKB integral as 

a/a 1 / 1 m r 



/-oo A _ - T 4 i/2 r ( J+n ^ V ' 

V 2+2n / 



Let us now consider the next order term in the 1/M expansion of (|2.12|) . There are 



two contributions to this constant piece. There is one contribution because ( |2.13| ) was 
integrated to oo instead of yo. Hence we should subtract from (|2.13| ) the term 



ly dV i^T = + 0(1/V/X) - 

The other contribution comes from integrating the integrand in ( |2.12| ) near yo. Subtracting 
off the leading order term, this contribution is given by 



A If" 1 ff'\ 2 l\ev 



yo 



oo 



/ 2 / 4 v /; v / 



f ^ [f^^ - l/x) - (l-I)(„ + l) + 0(1/VA). 

(2.15) 
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We have used the fact that /"// and /'// are almost constant near the turning point 
y = y for large VX. Hence, using O, (|2l3l) , (CT ) and (|2T5|) in (|2T2l) , we find that 

/ j " ( -+" ) \ 2 



16tt c 



T 2 m(m + n) + 0(m°) m > 1. 



(2.16) 



2+2n 



In the QCD3 case, n = 1 and hence the mass relation is 

M 2 = 87iYrf|^ T 2 m(m + 1) + 0(m°) m>l 
« 5.74216 (7rT) 2 m(m + 1). 



(2.17) 



We have factored out a 7r 2 term in the second line of ( |2.17|) to match the units used in . 
Table 1 compares the WKB expressions with the numerical results found in and we see 
that the agreement is very close.0 



m 


WKB 


Numerical 


1 


11.4843 


11.5877 


2 


34.453 


34.5270 


3 


68.906 


69.9750 


4 


114.853 


114.9104 


5 


172.265 


172.3312 


6 


241.171 


241.2366 


7 


321.561 


321.6265 


8 


413.436 


413.5009 



Table 1: Comparison of ++ glueball masses squared in units of 7r 2 T 2 . The WKB approximation 
is very close to the numerical results in 0], with a small difference approaching 0.064-7r 2 T 2 for 
large m. 

In the case of QCD^, we have n = 2 and hence the WKB relation 

r m 

M = 4tt 3/2 — ^-T v / m(m + 2) + 0(m _1 ) m > 1. (2.18) 

Table 2 compares the WKB results in ( |2.18| ) to the numerical results of [[T7||. Again, we 
find that the WKB result quickly approaches the numerical eigenvalues. 

1 In H a WKB expression was given with a numerical factor of 6. 
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m 


WKB 


Numerical 


1 


9.39 


9.85 


2 


15.3 


15.6 


3 


21.0 


21.2 


4 


26.5 


26.7 


5 


32.1 


32.2 


6 


37.6 


37.7 



Table 2: Comparison of glueball masses for QCD4 in units of T. The WKB approximation 
should approach the numerical result as a function of 1/m. 



We conclude this section by examining the behavior of the mass spectrum in the limit 
n — > 00. Taking this limit we approach p = 5, corresponding to six dimensional euclidean 
QCD with a compactified direction. In the large n limit, the mass equation in ( |2.16| ) 
reduces to 

M 2 = 4tt 4 T 2 (c + — ) m > 1, (2.19) 

where C is a constant to be determined. Thus, it appears that for nonzero C there is a 
finite gap in the spectrum, but above this gap the spectrum is continuous. We can see 
this more clearly by taking n — > 00 limit in (|2.6| ). In this case (|2.6| ) reduces to Legendre's 
equation and so the solutions that are regular at x = 1 are Pi(2x — 1) where A = —£(£ + 1). 
If A < 1/4 then Pi{2x — 1) is not normalizable at infinity. If A > 1/4, then Pg(2x — 1) 
is plane wave normalizable. Therefore, we find that the constant in ( |2.19| ) is C = n~ 2 
and thus there is a gap. This unusual behavior for the six dimensional theory is probably 
related to its nonlocal nature |T5],[TI|. 



3. Glueball masses for generalized supergravity backgrounds 

In the previous section we have seen that the analytic WKB expressions give accurate 
results for the eigenvalues of the dilaton wave equation. This strengthens are confidence 
in the procedure, and encourages us to use it in more general situations. 

In this section we discuss the WKB approximation for ++ glueballs in more general 
supergravity backgrounds. We will argue that that there is universality in the spectra 
which only depends on the dimensionality of the AdS space at infinity and the singularity 
structure at the horizon. Using results derived here we can find WKB approximations for 
^ glueballs in QCD4 and glueballs in QCD 3 . We can also explain the stability 
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of the ++ spectrum for models coming from rotating branes [13], and the change in the 
spectrum for the ^ glueballs in these same models. 

The only assumptions that we make are that there are angular independent solutions 
to the dilaton wave equation in ( |2.4| ) and that for large U the metric approaches an AdS 
solution. In this case, we can reduce Q2.4Q to 

du(f(U)du(l>) + M 2 h(U) = 0, (3.1) 

where 

f(U) = ^e- 2 *g uu , h(U) = Jge-^g™. (3.2) 

Let us suppose that there is a XJ$ such that near U = Uq, 

f(U) ~ (U-Uo) s h(U) ~ (U-U ) q . (3.3) 

The assumption that asymptotically the solution is AdS implies that f(U) ~ U 8 ~ p and 
h(U) ~UifU» U . 

Let us define e y = U - U , f(y) = f(e y + U )e sy and h(y) = h(e y + U )e qy . If we let 
-tp, then (pJ]) becomes 
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tfl> + V(y)1> = 0, (3.4) 
where 

v(y) = mv— m + i _ 1 3 («"-"'/(■>>), 

/(») 4 \ *-Vvf{y) J 2 c (-i)»/( y ) 
For large negative and positive y we have 

V(y) « CiM 2 e (9+2 - s)y - ^(a - l) 2 y«0 

4 (3.6) 

where C\ and C2 are unimportant constants for this discussion. Hence, we find two turning 
points which for large enough M can be approximated as 

1 , / 4CiM 2 

yi ~ ; ^ lo g 



q + s-2 °\(s-l) 



1 , / 4C 2 M 
2/2 ~ log 



( 3 - 7 ) 



5 — p \ (7 — p) 
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The WKB approximation is then 



1 



112 



m + -jir = J dy^V{y). (3.8) 



The leading order contribution to the integral in (|3.8| ) is 



Ml = M l + °°dyeyJ^M = M f°°dU^[^, (3.9) 



oc 



e s yf(y) Ju V 9 



where we have used (|3.2|) . It is clear that ( |3.9|) sets the inverse mass scale for the glueballs. 



We can compare this to the scale coming from the heavy quark potentials. In this latter 
case, the string tension along the brane, as a function of the energy scale U is 

<U) =^-g xx (U) (3.10) 

and at large quark separation the string tension approaches -^gxxiUo). Hence we can 
rewrite (|3.9|) as 

poo 

Mi = M \ dU^Jd{2ixo{U))- 1 ' 2 . (3.11) 

In other words, the inverse mass scale is a one-loop integral of the square root of the inverse 
tension integrated over all energy scales with a measure ^guu- 

Let us now consider the next to leading order corrections. The computations are 
similar to those in (|2.14 ) and ( |2.15| ). Using ( |3.6| ), the correction coming from the turning 
point at y = j/2 is 



(p-5)y/2 + 



POO 

- ^M 2 C 2 / dy e 
J" 2 dy (J M*C 2 elP-*)v - \{7 - p) 2 - VM^e^J = - (j^) |' 

while the correction coming from the turning point at y = y\ is 

- y^MKh / dye^ q+2 - s)y ' 2 + 

J — oo 

J^°° d V i^j M 2 C 1 eii+ 2 ~ a )y - i(s - l) 2 - ^/M 2 C 2 e( q+2 - s ^/ 2 j 



(3.12) 



s — 11 n 



q+2-s2 

(3.13) 



Putting everything together, we find that the WKB masses are 
- 2 < 2 Is -II 



7T 



^2 m [ m + 5 



+ 



+ 0(m°) 



m > 1. 



(3.14) 



p q + 2 — s 

By assuming that the supergravity solution is asymptotically AdS, we have chosen a 
particular singularity structure for the point at spatial infinity. The basic arguments used 
here are still applicable even if the solution is not AdS, so long as the singularity structure 
at infinity is known. 

We now consider some of the examples discussed in the recent literature. For the case 



of rotating nonextremal D4 branes considered in ||12| , p!3|1 , the dilaton wave equation reduces 
to 

d u [u (u 6 - (4gN) 2 a 4 u 2 - u%) d u <£\ + 4gNM 2 u 3 <p = 0, (3.15) 

where a parameterizes the angular momentum, and we have replaced U with u 2 = U, to 
match the form of the equation in [T3]. If a = this reduces to the nonrotating D4 brane 



equation in (|2.5| ). The horizon occurs at u = uq with 



Uq — a Uq — My 



0. 



(3.16) 



Hence, it is clear that (|3.15| ) has the form of (|3.1|) and (|3]3|), with s = l,q = 0, for all 
values of a. Therefore, we find that the WKB expression for the masses is 

-2 



m(m + 2) 



7T" 



4gN 



du 



u 



Uq 



^u 6 - (4gN) 2 a 4 u 2 - 



Un 



m > 1. 



(3.17) 



Since the singularity structure does not change when a is varied, we see that the WKB 
mass ratios do not change either. Thus, we see the reason for the stability of the glueball 



masses observed in ||13|| . This might also explain why supergravity glueball results are 
reasonably close to lattice results. 

We can compute the integral in ( |3.17| ) exactly in the two limits a = 0, oo. The result 
for a = is in (|2.18| ). For large a we find 

4 



n 



a 2 m(m + 2) 



m > 1. 



(3.18) 



The next examples are the ^ glueball masses for the rotating nonextremal D4 
branes. In this case, the equation of motion for one of the angular components of the R-R 
1-form field is f|13| 



d u [«V - (4 9 N) 2 a 4 )d uX \ +4gNM 2 J^J^^.e X 



u 6 - (4gN) 2 a 4 u 



0. 



Un 



(3.19) 
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This equation also has the same form as (|3.1|) , but with s = and q 
Therefore, using ( |3.14| )we find 

.2 



■1 for generic a. 



to(to + 3) 



7T 



4gN 

for the ^ masses. In the a 

M = 



du 



u 



V^ 6 - (4gN) 2 a 4 u 2 
= limit this reduces to 

r (2) 

4tt 3/2 — ^f|T v / m(m + 3) 

" V 6 ) 



6 



TO > 1, 



Ur, 



TO > 1. 



(3.20) 



(3.21) 



Table 3 shows a comparison of the WKB masses to the numerical results in [17] . The WKB 

result for the mass ratio between the lowest level and ++ states is M |_/M_|__|_ = 

2/v^3 ~ 1.155 which is reasonably close to the numerical result. In fact the difference 
between the WKB and numerical results is smaller than present day lattice errors. 



m 


WKB 


Numerical 


1 


10.8 


11.8 


2 


17.1 


17.8 


3 


23.0 


23.5 


4 


28.7 


29.1 


5 


34.3 


34.6 


6 


39.8 


40.1 



Table 3: Comparison of + glueball masses for QCD^ in units of T. 



However, in the limit that a — > oo, the WKB structure will change for the + states 



This is because the singularity structure of (|3.19|) changes. In fact, in the large a limit 
we end up with the same equation as the dilaton. Hence we find the same WKB masses. 
The only difference is that we have to discard the lowest eigenvalue [13]. Hence the WKB 
mass ratio in the large a limit is m++ = \/8/3 ~ 1.63. This is again close to the numerical 
result of 1.59 and it is also close to the lattice result of 1.61 ± .19. 

Our final example is the WKB spectrum for the 0~~ glueballs in QCD3. After a 
rescaling of the NS-NS 2 form field, the relevant component satisfies the wave equation 

du [U 5 (U 4 - Ufydux] + gNM 2 U 5 X = 0. (3.22) 
This does not have quite the same form as ( |0| ) because of the extra U 4 term, but the 
WKB analysis is almost identical and results in the spectrum 



M = VSn 



Ta/to(to + 3) to > 1. 



(3.23) 



Table 4 is a comparison of the WKB and numerical results []TIJ. Again, we find close 
agreement. Finally, the WKB mass ratio for the lightest states with 0~~ and ++ quantum 



numbers is 



M_ 
M. 



++ 



V2. 
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m 


WKB 


Numerical 


1 


4.79 


5.11 


2 


7.58 


7.82 


3 


10.17 


10.36 


4 


12.68 


12.84 


5 


15.16 


15.29 


6 


17.61 


17.73 


7 


20.05 


20.15 


8 


22.48 


22.57 



Table 4: Comparison of glueball masses in units of tyT. 



4. The type nonsupersymmetric model. 

In this section we make some general statements about the type model [ pO] , pT 



recently discussed in [jig] . We find asymptotic solutions for the lowest order supergravity 
approximation. We find a running of the effective coupling, but no confining behavior and 
no normalizable glueball solutions. 

The type model has a closed string tachyon, no fermions and a doubled set of R-R 
fields. In particular, there is no longer a self dual constraint on the 5-form field strength. 
Since the number of R-R fields is doubled, so are the number of D brane types [p2[ . Hence 
one can have D3 branes that are electric instead of dyonic. If we have N parallel electric 
D3 branes, then the low energy effective action on the branes is thought to be SU (N) 
QCD with adjoint scalar fields, but no fermions. Hence, there is no supersymmetry and 
the coupling will run. There is no open string tachyon [p2 |, so there is no tachyon in this 
QCD model. 



The authors in |T5[ argued that the closed string tachyon can get an expectation 
value, and that its mass squared gets a positive shift from the background 5 form flux. 
The background tachyon field acts as a source for the dilaton, so the dilaton is no longer 



constant. One then makes the following ansatz for the metric |L5[ 

ds 2 = e** [e^'^dp 2 + e -Udxl + e^dflfj , (4.1) 
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where <f>, £ and r\ are functions of p only. The equations of motion then reduce to a Toda 
like system with an action [p!5|B 



S = dp 



V(<f>,£,ri,T) = -T 2 e^ + ^- 5 ^ + 20e~^ -Q 2 f- 1 (T)e~ 2 ^ 
2 



(4.2) 



and a constraint 

\fi + \? + \T 2 - 5f + C v) = 0. (4.3) 

Q is the total D3 brane charge which is proportional to N, T is the tachyon field and f(T) 
is a function given by fT5 



f(T) = 1 + T+-T 2 + 0(T 3 ). (4.4) 

For large Q, the tachyon expectation value is determined by setting f'(T) = 0. As a first 
approximation, we may assume that the tachyon is constant as a function of p. 

If T = then the solution reduces to the M — 4 solution. When T is nonzero, then 
all fields are coupled and there is no known analytic solution. However, we can attempt to 
find approximate solutions that are valid in the UV and IR regions. In the UV, we expect 
the dilaton field to be relatively constant, at least compared with £ and r\. Assuming that 
is constant and thus ignoring its kinetic term, the equations for £ and r\ can be solved 
exactly, at least in the near horizon limit. In this case we find 

e ? = C x p e" = C 2 p 1/2 

1 n i A C) /2 2Q 2 

4 C| + C'if(T) 1 " ° (4.5) 

Ve**^ + ^-5 = 

One can easily check that this satisfies the constraint equation in 04.3j ). If we plug this 
back into the metric, we find that the solution is still AdS$ x S$, but the curvatures of 
the two spaces no longer match, £5 now has smaller curvature then AdS§. In this case the 
Ricci scalar for the total space is proportional to 

R ~ T 2 eK (4.6) 



We are using units where a' = 1. 
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Using the £ and r\ solutions as inputs, we can go back and find an approximate solution 
for 4> m terms of p. Using the ansatz = Co(log(p/ po)) a , and plugging this into the 
equation of motion for <fi 

+ br 2 e^+^-5v = q 5 (4.7) 

we find that the ansatz is a leading order solution if a = — 1 and Cq = — 8C|/(T 2 -\/C i i)- 
po is an integration constant and we assume that po >> 1 in order that the gauge theory 
length scale is much greater than the string scale. Setting p = w -4 , and using the lowest 
order solutions for C\ and Ci from ( |4.5| ), we learn that the leading order behavior for the 
coupling is 

i QT 4 ( log ±) 

e -* = _L = - V (4 8) 

9ym 4096v^J(T) 

Thus we find a running coupling, but instead of a linear log dependence, the coupling runs 
with a log squared! One can easily check that to leading order in 1/ log u, the constraint 
equation is still satisfied. We can also estimate the range of validity for this solution. 
Computing the leading order corrections to C\ and C%, one finds that 

We can also compare the terms in the potential that depend on the tachyon. Since 

„ Q* r \T)e-* ~ ^, (4.10) 

our solution with constant T and f'(T) = is valid so long as log(u/wo) >> 1- 
The metric in the large u limit is 

Hence we can trust the supergravity solution only if T « 1, since log(u/wo) >> 1- 
However, it is clear that T ~ 1 if f'(T) = 0, hence one should not expect the supergravity 
result to be particularly trustworthy. Indeed, we have found that while the supergravity 
computation results in a coupling running to zero, it runs with the wrong power of log it. 

Nevertheless, the effective coupling between a heavy quark and its antiquark does 
appear with the expected log dependence. Using the Wilson line computation of p3| , p4] , 
one finds that the quark potential is given by 

256 7r 3 
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where we have plugged R 2 = T2 i ^ u / UQ ^ into the expressions derived in p3| , |24]1 . L is some 



length that can be adjusted to be much longer than the string scale. Recall that the M = 4 
potential comes with a coefficient y/gN. It is this coefficient, and not the 't Hooft-Polyakov 
tension gN, that plays the role of 1/a' for the supergravity models. For a string theory 
with extrinsic curvature, the string tension has a (log it) -1 dependence [PSLEB]. Hence our 



result is in line with rigid string results, so long as one remembers that the tension is y/gN. 

Even though we cannot really trust the supergravity solution for large u, we might 
be able to trust it for small u. However, here we will see that the situation is even worse. 
In particular, we will find that the dilaton wave equation in the lowest order supergravity 
background has no normalizable glueball solutions. Moreover, the heavy quark potential 
is not confining. 



In order to study the IR behavior, let us follow the suggestion of |T5j and search for 
solutions for the to da system assuming that the second term in V is small compared to 
the other terms. This corresponds to a small curvature for the S$. Dropping this second 
term in V and assuming a constant T, one can now find an exact solution to the equations 
of motion that satisfies the constraint. The solution is 

= C 2 p 5 / 9 e*> = C 2 p 5/9 e* = jg_ p, (4.13) 

with the relation 

20(2/(T)) 1 / 4 C| -9T 2 y/3QC Q = 0. (4.14) 

Comparing all terms in V, one has e^ + ^ _5?7 ~ p -2 , e~ 2 ^ ~ p -2 , but e~ 4r] ~ p -20 / 9 . 
Hence this solution is valid for large p. From (|4.13| ), the coupling blows up as p — > oo and 
after substituting p = 1/w 4 the metric is 

ds> = i» (l6^ + ^MSm^i + du-VW^j . (4.15) 

C*2 remains as a leftover integration constant and is ultimately determined by matching to 
the UV solution. 

In the small u limit, the T dependent terms in the potential are now comparable, so 
the tachyon expectation value is no longer at f'(T) = 0. Instead, plugging in the solution 
of (|4.13| ) and Q4.14j ) into the tachyon equation of motion, one finds that T = if 

10/(T) + Tf'{T) = 0. (4.16) 
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Using the function in (|4.4| ) for f(T), one finds that there are no real solutions of (|4 . 1 6 ) for 
T. However, if we were to include the cubic term in f(T) or include the quartic term ciT 4 
in the action, then a real solution would exist. 

As in the UV, the curvature in the IR is small if T « 1. However, we expect solutions 
to (|4.16|) to be T ~ 1. The curvature is now at the string scale, so we cannot truely trust 



the supergravity solution in this limit either. 

We can attempt to find glueball solutions. Using the asymptotic expressions in the 
dilaton equation of motion we find the equations 

2Q M 2 

d u [u 5 d u (f\ + j== u(f) = u » u 

2 (4.17) 

du [u5du(p] + ^wm u21/9<p = u<<uo 

For large u, we see that the dilaton equation of motion reduces to the Af = 4 result. It 
appears that the only other singularity is at u = 0. Using the arguments of the previous 
section, one learns from ( fl.l7| ) that this singularity has s — q = 5 — 21/9 > 2, hence no 
glueball solutions are possible. 

We can also easily see that the potential between the heavy quarks does not confine 
for these solutions. If we define a new variable v such that 



9 

then the metric in ( f4.15|) is 



{ 3Q J 



1/2 



u A '\ (4.18) 



,2 180 / dv 2 , a 2 ; (C 2 /3) lo y27(r) ^2 , 
ds = T2 \ -2 + v dx \\ + g2-2 I • 1 1-1-9) 

From this metric, we see that R 2 — ^r, and so the heavy quark potential is 

1440 7T 3 

V « - 7^7 \ L » L Q . (4.20) 

r(|) 4 T 2 L v ; 

Hence, the supergravity result implies that the effective heavy quark coupling increases 
when going from the UV to the IR, but the quark potential does not develop a linear term 
and remains coulombic. 
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We could certainly generate a linear quark potential by going to finite temperature 
since the supergravity background would now have a horizon at finite u. The down side of 
this scenario is that the the theory will essentially be reduced to QCD in three dimensions. 

In the end one probably has to consider the full cr-model in order to get a linear quark 
potential and normalizable glueball solutions. At the very least, one could try including 
the a' 3 R 4 terms in the action to see if this would qualitatively change the behavior in 
the infrared. Perhaps one could then combine the UV results found here with IR results 
derived from the cr-model to say something concrete about the WKB glueball masses. 

5. Conclusions 

We have seen from the WKB mass expressions that there is some degree of universality 
for glueball mass ratios between different supergravity models. In the examples where the 
ratios change, the behavior can be attributed to a change of the singularity structure at 
the horizon. 

The finite temperature models do not exhibit a running of the coupling in the UV. 
However, we have shown that the type model has the desired behavior. Unfortunately, 
it does not appear to be confining in the IR. Hopefully confinement will appear in the 
solution for the full a-model. Or perhaps a model can be found that combines the desired 
features of the nonsupersymmetric models discussed here. 
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